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Abstract 

We have formulated the quasipotential method for the calculation of the rela- 
tivistic and radiative corrections to the magnetic moment of the two-particle bound 
state in the case of particles with arbitrary spin. It is shown that the g- factors of 
bound particles contain 0(a 2 ) terms depending on the particle spin. Numerical 
values for the g-factors of the electron in the hydrogen atom and deuterium are 
obtained. 
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The study of electromagnetic properties of hydrogen-like atoms and ions in quantum 
electrodynamics (QED) is one of the basic tasks in the theory of two-particle bound 
states. The experimental verification of the g-factors calculation for the bound particles 
is carried out during many years [[I], EL The values of the electron g-factors in hydrogen 
atom, deuterium, tritium and helium ( 4 ife + ), measured in many experiments, are in good 
agreement with theoretical results. Recently the scope of the experimental investigation 
of hydrogenic ions was essentially inlarged 0, |J. These experiments generate a need for 
the new theoretical study of different contributions to the gyromagnetic factors of the 
bound particles || [pj, [Tj. At present, the most accurate value for the electron g-factor is 
obtained in the experiment with the hydrogenic carbon ion 12 C 5+ (Z=6) || ||: 

g e e xp { 12 C 5+ ) = 2.001 041 596 4(8)(6)(40), (1) 

where the statistical (8), systematical (6) errors and inaccuracy, connected with the elec- 
tron mass (40), are shown in brackets. Theoretical investigations of the electromagnetic 
properties of the hydrogen- like atoms, which were done in [8-14], showed that the gyro- 
magnetic factors of the bound particles can be written in the from: 

g{H - atom) = 2 + Ag rel + Ag rad + Ag rec + ... (2) 

Relativistic corrections Ag re i, radiative corrections Ag ra d, recoil corrections Ag rec in (2) 
were calculated with the accuracy up to terms of order a 3 (m/M) and a 2 {m/M) 2 in ||, |9] 
on the basis of the quasipotential method for spin 1/2 particles, composing the bound 
system. The dots designate other possible contributions to the g-factor. In addition, due 
to the experiments with deuterium, hydrogen-like ions, which have the nucleus of arbitrary 
spin, there is need to extend the calculational methods of the g-factors on this case. In 



paper 15 it was suggested an approach for the calculation of the corrections to the 



gyromagnetic factors, based on the Bargmann-Michel-Telegdi (BMT) equation [IS]. The 



conclusion about independence of the binding corrections on the magnitude of the spin 
of the constituents was also formulated here. In this work we extended the quasipotential 
method for the calculation of the magnetic moment of two-particle bound state to the 
case of arbitrary spin particles and calculated main contributions to the corrections in 
eq.(2). 

The interaction of the massive particles of arbitrary spin with the electromagnetic 
field is studied on the basis of different methods during a long time [17-25], but still this 



problem is far from its final solution. It was shown in [18] that the particle of arbitrary 
spin must have at tree level approximation the gyromagnetic factor g=2. In general case, 
the matrix element of the electromagnetic current for the particle of arbitrary spin s is 
determined by means of (2s+l) form factors (charge, magnetic, quadruple, et al.). When 
studied the magnetic moment of the simple atomic systems, it may be possible to take 
into account the form factors of the minimal multipolarity, describing the distributions 
of the electric charge and magnetic moment. The one-particle matrix element of the 
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electromagnetic current operator between states with momenta p and q can be written 
as follows: 

J, = U(p) {iVFf + J-£^FF 2 P } U(g) (3) 

The wave function U(p) of particle with arbitrary spin, entering in (3), can be presented 
in the form (see, for instance 124], 125!): 



U 



v 1 Ia\6i2...a v , 



p + q = 2S, 



(4) 



where spin-tensors £, 77 are symmetrical in upper and lower indexes. For the particle of 
half-integer spin p=S+l/2, q=S-l/2. In the case of integer spin p=q=S. The Lorentz 



transformation of the spinors £ and r\ can be written in the form p5| ^6 

£ = e 2 £ , V = e 2 £0, 



(5) 



where the direction of the vector coincides with the velocity of the particle, th<f> = v. 
The generator of the Lorentz transformation E is equal: 



P+Q 



i=l i=p+l 

and <7j acts on the ith index of the spinor £0 as follows: 



(6) 



(7) 



In the standard representation, which is introduced in analogy with the spin 1/2, the free 
particle wave function (4) may be written with the accuracy (f/c) 2 in the form: 



U(p) 



1 + 



mi 

Sm 2 
Eft 



£0 



/ 



The components of the matrix E MV in (3) are the generators of the boosts and rotations 



2§: 



J n0 



£„ 
-E„ 



-2ze 



rank 



Sk 
s fc 



1 2S 

s= ~£°v 
z i=i 



(9) 



The general expression for the magnetic moment of the two-particle bound system reads 

as a m-. 



M — —- 



I 



d_ 
dA 



x <K A \J{0)\K B > 



A = K A - K 



(10) 
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where the matrix element of the current operator between the bound states 

<K A \J,(0)\K B >= J ^^5(p 1 +p 2 -K A )^ A (p)T,(p,q,E A ,E B )x (11) 



dq\dq 
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k b w»w ' «« (27r)3 > 

is expressed by means of the wave function of the bound state \I> ^ b (p) and the generalized 
vertex function T M presented in Figure 1. The vertex function is determined through 
the five point function: 

R p =< 0\Mt,xi)MtiX2)JMMr,yi)MT,y 2 )\0 >, (12) 

projected onto the positive-energy states, 

r„ = g-^mct 1 , 4+) = utu*R„u 1 u 2 , (13) 

where G is the two-particle Green function. We study the loosely bound composite 
system, so it is possible to make the perturbation expansion of all quantities T, R and 
G~ x in particle interaction. 

r = r(°) + r« + ..., r = r + r 1 + ..., g~ x = gz 1 - v x - ... , (14) 



r(°) = G^R G \ (15) 

r« = Gv 1 r 1 Gv 1 -v 1 g o tM-tWg v 1 ,..., (ie) 

where G is the Green function of two noninteracting particles and V\ is the one-photon 
exchange quasipotential (see eq. (19)). 

The transformation law for the wave function \1/ g g (p) of the system of bound particles 
with spins S\, s 2 from the rest frame to the reference frame, moving with the momentum 
K B , was obtained in |27| : 



5(p 1+ p 2 - K b )^ Rb ($) = Df\R w )D^(R w )Jj-^%{^)5ipl+p° 2 ), (17) 

where D S (R) is the well-known rotation matrix and Ryy is the Wigner rotation associated 
with the Lorentz transformation A^ b . (E, K B ) = A^ s (M, 0); (e,p) = A(e°,p i ), E = 



M 2 + K B , e(p) = y/p 2 + m 2 . The exact expression for this rotation matrix is the 
following ||: 

D S (R W ) = S-'ifiSiK^Sip*), (18) 
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(b) 

Figure 1: Generalized two-particle vertex function T^: Diagram (a) represents Di- 
agram (b) describes where the bold line denotes the negative-energy part of the 
propagator. 
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where S(p) is the Lorentz transformation matrix of the spinor wave function (4). The 
quasipotential bound state wave function ^o(p°) in the rest frame of the composite system 
satisfies to the quasipotential equation [[28] : 

where the relativistic reduced mass 

E X E 2 _ M 4 - (m\ - mf) 2 _ M 2 - m\ x + m\ 2 

m ~ ~W ~ 4AP ' 1,2 ~ 2M 

M = Ei + E 2 is the bound state mass, 

[M 2 - (mi + m 2 ) 2 ][M 2 - (mi - m 2 ) 2 ] 



6 2 (M) 



AM 2 



In the nonrelativistic limit this equation reduces to the Shroedinger equation with the 
Coulomb potential. The D s functions in (18) can be obtained in the approximate form, 
using (5): 

s gMggxgg) k 2 -(M b )(M b ) p>k b -{2f){2k b ) 

D ( R ^)- 1 + 4^ + + • (20) 

The main contribution to the vertex function r M , which is determined by Feynman dia- 
gram (a), shown in Figure 1, can be expressed as follows (in the Breit frame): 

f (°>(£ q) = ^(pi)eijf ! + [5i x a] }[/i(gi)5(p 2 - <f 2 ) + (1 <- 2), (21) 
where Ef 2 (0) = ei >2 , -Fi P 2 (0) = ei j2 /ti j2 , and the matrix 

is the natural generalization of the corresponding expression of the Dirac matrices 7 for 
the spin 1/2 in the standard representation. To simplify different terms of eq. (21) it is 
helpful to use the following relationships 



[£i, Sj] = 4ie ijk s k , [£ i; Sj] = ie ijk E k . (23) 

To construct vertex function T^°'(p, q), accounting for the terms (v/c) 2 , we use explicit 
expression of the wave function (8), transforming the different parts of the matrix element 
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(21) by means of the equation of motion for the spinor U(p). The following relations are 
valid (taking into account 5(p 2 — q 2 )): 



m + oi e 2 - ip° (Si[p° x A] 

Ui{pi) V -±P±U{q x ) = 2p° - -|A + V 2 

2mi M m{ 



(24) 



- ,^ £i(pi) -ei(gj) -p 



' X p°], 



m 



(25) 



?7 



iCpi)^! x =gxA]-ij {p° (^[p" 8 x A]) + [p° x S^A)} . (26) 



Bound state effects in the vertex function r M are determined by the diagram (b) of Figure 
1. Accounting for iteration terms with the quasipotential, we can represent the corre- 
sponding expression in the form B M: 



f^(p,q) = U^p 1 )U;(p 2 )^\M{-\^)B 1 B 2 V(q 2 -p 2 ) + 



(27) 



+B l B 2 V{p 2 - g 2 )A^Vi)^l )Ui{qi)U 2 {q 2 ) + (1^2) 



V(k) = B 1 B 2 { ( 1 + ^-f t k ) ( 1 - -P-T 2 k ) 



2m! 



2m 9 



:28i 



A\ + ^fASi x k]) (A 2 - ^-B 2 i[S 2 x £]) J^, 

where the negative-energy projection operator A~(p) « (1 — B)/2 — Ap/2m, Ki )2 are 
the anomalous magnetic moments of the particles. The matrices Ai )2 , B% y2 are also the 
generalizations for a\ }2 , /3i j2 used in the case of spin 1/2 particles as in (29): 



A 



E 
E 



, B 



I 
-/ 



(29) 



Both parts of the quasipotential (28) give contribution to the magnetic moment of the 
system with the accuracy (v/c) 2 . Substituting (27), (28) and (8) in (11) and calculating 
derivative in A to eq. (10), we obtain: 



M 



dp^(p)-ii_/2(l + Kl )5i[l + iVi + N 2 ] + (1 + W PX ^ 2 XP1] + (30) 



+ (1 + K 2 



(2tt) 3 7 "'" uvr, 2ei 

[p x [s 2 x p]] Ef e 2 (p) 



2m? 



m\m 2 



M 



1 + Ax + A 2 + 



(M - ei - e 2 ) E? 



m 2 





\ d' 








op 



+ 
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where 



Ni 



2m? 



In the case of S-states the expression (30) may be essentially simplified: 

-* 1 ei _ 1 e 2 

•M = -Pi bound < Si > +-<?2 bound < «2 >, 

2 mi 2 77i2 



where the g-factors of the bound particles are equal: 

3/ci 



#1 bound — #1 % 1 



3m? 



2(1 + Ki) 



+ 



< p 2 > 
2m? 



i-<S> + 



(31) 



(32) 



(33) 



+ 



mf 
m 2 . 



< £ 2 > 



+ 



e 2 < F 2 >< £ 2 > 
ei 3m 2 3 



< p 2 > 



[1 + K, 1 )6m 1 (m 1 + m 2 ) 



1 - 



e 2 mi 
ei m 2 



#2 feound — #1 bound (1 ^2), 



2^1,2 = 1 + «1,2- 



For the hydrogen-like ion (1 is the electron, 2 is the nucleus) we have: e\ = — e, e 2 = Ze, 

< p 2 >= mlml(Za) 2 / (m 1 + m 2 ) 2 , 



K S1 = 



<£?> 



1, K a 



<£ 2 > 



^p, s 2 is the integer nucleus spin 
4s2 „ +1 , s 2 is the half — integer nucleus spin' 



so the g-factors of the bound electron and nucleus take the form: 



Qe bound Qe \ 1 



2m|' ' 



S2 



m 2 (Za) 2 
3(mi + m 2 ) 2 



3^1 



2(1 + «i) 2 
mi 



2(mi + m 2 )(l + Ki) 



1 + 

m 2 J 



9N bound — gN\ 1 



m?(Za) 2 
3(m x + m 2 ) ; 



3K2 -gq-K ^ 

2(1 + « 2 ) 2 U As2 ' 



3 ml 
"2m? 



si 



3Z 



K S1 ) + 



m 2 



1 + 



m 2 



(34) 
(35) 



(36) 



2(mi + m 2 )(l + k 2 ) V Zmi, 

The expressions (35)-(36) were obtained from interaction amplitudes, shown in Figure 1. 
They contain corrections of order 0(a 2 ), 0(a 3 ), connected with the bound state effects. 
As it follows from (30), some of these contributions can be considered as relativistic 
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corrections for spin s particles (terms proportional to iVj). Other corrections refer to the 
spin-orbit interaction. Our calculation of relations (35)-(36) shows, that the terms of 
order 0(a 2 ) in g e bound, 9n bound depend on the spin of second particle - nucleus. This 



conclusion differs from the results of the paper [T5\, where such dependence is absent. In 
the case of spin 1/2 particles the obtained expressions (35) and (36) for the bound state 
g- factors coincide with the results of || [TO], [TTj, £5(J. The electron g-factors in hydrogen 
atom, deuterium, tritium as well as their ratios are very important from the experimental 
point of view JTJ. Experimental value of the ratio g e n/de d, obtained in JTSJ with high 
accuracy 

r 9e 



„exp 



H 



9e D 



exp 



1 + 7.22(3) ■ 10" 



(37) 



The theoretical expression for this ratio can be written from (35) as follows: 

25 m 2 







th 


r th = 


9e H 


= l + a 2 




_9e D _ 






1 m\ 
16 m 2 . 



(3? 



As we pointed out above, the new approach to the calculation of the various order 
contributions to the magnetic moment of the loosely bound system was suggested in the 
work [[15]. This method is based on the relativistic semiclassical equation of motion for 
spin. Constructed in jl5] on the basis of this equation the interaction Hamiltonian of 
particles of arbitrary spin with the external electromagnetic field leads to the g-factors 
of the bound particles, which are independent of its spin. It is known that the BMT 
equation is approximate one: it is linear on the particle spin, field strength F uu , which 
doesn't contain coordinate dependence. When the spin s particle is in the bound state in 
the external homogeneous magnetic field, some terms, omitted in the approximation of 
the BMT equation, can give definite contribution to the g-factors of the bound particles. 
In this work during the calculation of nucleus spin- dependent terms in the gyromagnetic 
factors of the bound particles, composing hydrogen-like ion, we used the approach, pro- 
posed in the papers [|24], 155] , for the description of the electromagnetic interactions of 
particles with arbitrary spin. New contributions to (35), (36), (38), as compared with ||, 
were obtained after the replacement of the ordinary boost generators a in the case of spin 
1/2 particles by operators (29). The value of the correction in r th , connected with the 
deuteron spin is equal to Ar th = 5a 2 m 2 /72m 2 = 0.001 • 10~ 9 (1=1, Z=l, m nuc i = 2m2 (m.2 
is the proton mass). It lies within the limits of the experimental error, as it follows from 
eq. (37). Numerically, the ratio (38) r th = 1 + 7.237 • 10 -9 , which is in good agreement 
with (30). Nucleus spin dependent corrections in (35), (36) are the functions of Z and 
the number of nucleons in the nucleus N. Despite the increasing of these corrections with 
Z as Z 3 the growth of the nucleon number N in nucleus leads to opposite effect. So, in 
the case of ions with nucleous spin 1^0 the numerical value of these corrections is out of 
experimental accuracy reached for 12 C 5+ . For the ion 12 C 5+ 1=0, so Ki = and the corre- 
sponding spin-dependent correction is equal to zero. At present time the measurements of 
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the electron g-factors in the ions 16 7+ and 32 S 14+ were carried out [Q, but their nucleus 
have also spin 1=0. From our viewpoint, it will be interesting to measure the values of 
electron g-factors, using the Penning traps || |4[] for such ions, which have nucleous spin 
1^0 and the ratio Z 3 /N 2 would be reached large values. One of the such ions is the ion 
59 Co 26+ , which has I = 7/2, Z 3 /N 2 pa 5.65, and the value of spin-dependent correction in 
(35) is equal to 0.1 • 10~ 9 . As was pointed out by W. Quint [Q, the measurement of the 
bound electron g-factors with the accuracy exceeding 1 ppb may be realized in the near 
future. 

We are grateful to W. Quint and V.M. Shabaev for the information about new experi- 
mental data on the g-factors of bound particles, and to S.G. Karshenboim, LB. Khriplovich 
for useful discussion of the questions, connected with the magnetic moment of the two- 
particle bound state. The work was performed under the financial support of the Russian 
Foundation for Fundamental Research (grant 00-02-17771) and the Program "Universities 
of Russia - Fundamental Researches" (grant 990192). 
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